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ABSTRACT 

We use extensive multi-wavelength photometric data from the Great Observatories Origins Deep 
. . . Survey to estimate photometric redshifts for a sample of 434 galaxies with spectroscopic redshifts in 

CO I the Chandra Deep Field South. Using the Bayesian method (Bem'tez 2000), which incorporates red- 

■ shift/magnitude priors, we estimate photometric redshifts for galaxies in the range 18 < Rab < 25.5, 

I giving an rms scatter cr( ^^^^ — ) < 0.11. The outlier fraction is < 10%, with the outlier-clipped rms 

being 0.047. We examine the accuracy of photometric redshifts for several, special sub-classes of ob- 
jects. The results for extremely red objects are more accurate than those for the sample as a whole, with 
Jjfi i a = 0.051 and very few outliers (3%). Photometric redshifts for active galaxies, identified from their 

X-ray emission, have a dispersion of tr = 0.104, with 10% outlier fraction, similar to that for normal 

04 , galaxies. Employing a redshift/magnitude prior in this process seems to be crucial in improving the 

agreement between photometric and spectroscopic redshifts. 

^ . Subject headings: galaxies: evolution — galaxies: spectroscopy — galaxies:photometry 

00 ; 

^ ■ 1. INTRODUCTION in GOODS are fainter than the practical hmit for spec- 

5 ' The successful installation of the Advanced Camera troscopy. ^ Accurate photometric redshifts are therefore 

Or c ^Ar^c^ u j j-u tr uui c i essential, m order to estimate distances to galaxies host- 
tor Surveys (ACS) on-board the Hubble Space ielescope . , t . . i Ti 
rr^ /Tjorpx J A • r u J- • ri 4.1, c mg type la supernovae, to convert observable properties 
"i (HSl ), and the promise of new observatories like the Space / . , ^ { . .^ • ,, . r ^ 

O- T c J rp 1 -i-i /cTT3rrT7\ u J (sizc, color, magnitude) to those m the rest-irame, and 

^ Infrared telescope I'acility (SlRil , have opened new op- ^ ' ' ° ' , „ , , ,. ' , . 

' i -i- r u- 1 -f;! 1- ; - • to cxplore mauy other aspects of galaxy evolution. Ihis 

, portumties tor multiwaveiength surveys 01 the distant um- , ,, -11 r ^ ^ 1 . 

TT ifiu i-jj-j-j-j-T, letter provides details 01 the procedure used to estimate 

~ 1' verse. However, most 01 the galaxies detected m these ^ , ^ alaxi in the GOODS CDF S nd t 

' ' surveys are too faint for spectroscopic observations. The ''^ ^ /, ^ *^ ga axies m ^ , an o as- 

O ' u 4. 4. • J I,- -(-4- 4. I, • -4.1, 1 4.- 1 i sess the accuracy with which these are measured, using 
iZj • photometric redshiit technique IS the only practical way to , „ / . , . , , -n 

M , -4-1 1 • J J r J. 4.- 4.- 1 a sample 01 galaxies with spectroscopic redshitts. We 

estimate distances to these galaxies, needed tor statistical , , , ,, r 1 , • 1 1 -r- r • 1 

c/j ■ , J. r . 1 • . u • u 4-1 J also test the accuracy of photometric redshifts lor special 

rri , studies ot their properties. Ihe technique has the advan- , , ri-/,i.i .-i ■ .-r- ■ . 

, r .J. J 1 r 1 1 r r • 4- 1 subsampics ot obiccts that havc particular scientifac luter- 

• • ■ tage ot providing redshitts tor large samples ot taint galax- , .,, . ^v^^t^o , , : , , 1 1 • , 

>- , . .,1 w 1 J 4. • i 4. r 1 • 4-- est withm the GOODS data set: extremely red obiects 

les with a relatively modest investment ot observing time, fERO ) nd tiv al ti n lei fAGNs) 

but has the disadvantage of having coarse wavelength res- ^ ^' ^'^ ^^'^ ^ 

olution 1000 A) compared to a much higher resolution 2. the sample 

5^ , obtained spectroscopically 10 A). 

The Great Observatories Origins Deep Survey 
(GOODS) is designed to reach unprecedented depths in We limit our present analysis to the GOODS CDF-S 

two fields, the Chandra Deep Field South (CDF-S) and field, where the most complete, multiwaveiength photo- 

Hubble Deep Field North (HDF-N), using HST, Chan- metric data are currently available. The imaging data 

dra, and SIRTF. The majority of the objects detected and photometric catalogs are described in Giavahsco et 

al. (2003). Here, we have used ground-based optical 
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{U'UBVRI) and near-infrared {JHK,) data from ESO 
facilities (2.2m-WFI, VLT-FORSl, NTT-SOFI, VLT- 
ISAAC), as well as from HST/ACS (-B435V'606«775-285o)- 
Most of the imaging data cover the whole GOODS field 
of view, but a few (ISAAC JHKs, and FORSl BJ) do 
not. In total, there are as many as 18 independent pho- 
tometric measurements for each galaxy. We matched the 
point spread functions (PSFs) of the images, including the 
ACS data, as described in Giavalisco ct al. (2003). An R- 
band selected catalog was created with SExtractor (Bcrtin 
& Arnouts 1996). The photometry used for the photo- 
metric redshifts was measured through matched, 3 arcsec 
diameter apertures in all bands. 

2.2. Spectroscopic Observations 

The spectroscopic observations of CDF-S galaxies are 
from two separate samples: 

K20: A portion of this if-selected survey (^^(Vcga) < 
20; K,{AE) < 21.85; Cimatti et al 2002a, 2002b) lies 
within the CDF-S, and has reliable spectroscopic redshifts 
for 271 objects, with a median redshift (z) = 0.85. Using 
their spectra, the galaxies were classified as normal (i.e., 
emission and absorption line systems) or AGNs. 
GOODS/FORS2: The GOODS team is carrying out an 
extensive spectroscopic campaign in the CDF-S as an ESO 
Large Programme (C. Ccsarsky, PI). Here, we use data 
from the first observing season, obtained with the VLT- 
F0RS2 spectrograph using the 3001 grism {R « 860). The 
primary targets were faint {zs5o < 24.5), red (1775 — ^:850 > 
0.6) galaxies selected from the ACS imaging. A hetero- 
geneous sample of other galaxies was used to fill out the 
multislit masks. We use 163 redshifts that were measured 
with good confidence, with a median redshift {z) = 1.05. 

3. PHOTOMETRIC REDSHIFT TECHNIQUE 

We matched the spectroscopic and photometric cata- 
logs, and used the multicolor data to estimate photo- 
metric redshifts for galaxies in the spectroscopic sample. 
The comparison between photometric and spectroscopic 
redshifts then gives an estimate of the photometric red- 
shift accuracy. Because the two spectroscopic samples 
are selected differently, we have analyzed them separately 
and compare the results to explore any biases induced by 
color/magnitude selection criteria. 

We have performed extensive experiments with two vari- 
ant methods of photometric redshift estimation: spectral 
template minimization (Puschell et al. 1982), and the 
Bayesian method of Benitez (2000). In the technique, 
the photometric redshift of a galaxy is estimated by com- 
paring its multi-band photometry with spectral energy dis- 
tribution (SED) templates of galaxies with known types, 
shifted in redshift space, and integrated through the band- 
pass throughput functions. For each template, at each red- 
shift, the statistic is estimated, with the best-fitting 
redshift and spectral type found from the minimum 
value assigned to the galaxy. 

The Bayesian approach considers the distribution 
p{z\C\m), i.e. the redshift probability given not only the 
observed colors of a galaxy, C, but also its magnitude, m. 
This can be written as (Benitez 2000): 

p{z\C\m) (X p{z,T\m)p{C\z,T) 

BPZ is available from http://acs.pha.jhu.edu/~txitxo/, Benitez S 



The first term on the right is the redshift/magnitude 
prior, which contains the probability of a galaxy hav- 
ing redshift z and spectral type T, given its magnitude 
m, while the second term is the redshift/type likelihood, 
piC\z,T)^expi-xiz,T)y2.). 

Note that using a prior of this kind is not equivalent to 
assuming a particular luminosity function; the prior de- 
scribes the expected redshift distribution for galaxies of a 
certain magnitude, but does not include information about 
the galaxy magnitude distribution. The shape of the prior 
is estimated empirically from the HDF-N observations, as 
described in Benitez (2000). The best estimate of the red- 
shift, zb, is then defined as the maximum of p{z\C,m). 

We compared results from several different software im- 
plementations of the method, as well as different choices 
of the SED templates (empirical or synthetic), and found 
them to be in generally good agreement with one another 
and with the Bayesian Photometric Redshift (BPZ) results 
up to z ^ 1. However, at higher redshifts, the BPZ ap- 
proach gives significantly better performance, as described 
in §4, below. We have used the BPZ software which of- 
fers both simple minimization, as well as the Bayesian 
estimate. We use template libraries, consisting of E, Sbc, 
Scd and Im SEDs from Coleman, Wu & Weedman (1980) 
and two starburst SB2 and SB3 templates from Kinney et 
al. (1996). A two point interpolation between each pair of 
templates in the color-redshift space is performed, signif- 
icantly improving the redshift resolution. 

The BPZ uses a quality indicator, the Bayesian ODDS, 
which can be efficiently used to discard those objects with 
unreliable photo-z's. The ODDS is defined as the in- 
tegral of the probability distribution p{z\C',m) within a 
0.27(1 + zb) interval centered on zb- The ODDS would 
thus be 0.95 for a Gaussian P{z\C,m) with width a = 
0.067(1 -|- Zb), which is the empirically measured accuracy 
of BPZ in the HDF-N. 



4. RESULTS 

4.1. Overall Photometric Redshift Performance 

We quantify the reliability of the photometric redshifts 
by measuring the fractional error for each galaxy, A = 

hili£L-l£p£^ ^ W/g examine the median error, (A), the rms 

scatter ct(A), and the rate of "catastrophic" outliers, 77, de- 
fined as the fraction of the full sample that has |A| > 0.2. 

Table 1 presents the comparison between photometric 
and spectroscopic redshifts for both the K20 and F0RS2 
samples, using the complete photometric data set avail- 
able for each object, and for different combinations of 
magnitude limits and the ODDS parameter. Overall, the 
BPZ method yields o-(A) « 0.11 and 0.07 for the F0RS2 
and K20 samples, respectively. At fainter magnitudes 
{Rab > 25), the BPZ method gives smaller rms scatter 
than the conventional x^ minimization technique^ with no 
priors. The BPZ results are presented in Figure 1. The 
agreement between the photometric and spectroscopic red- 
shifts is excellent, with median offset (A) = —0.01. There 
is no substantial trend in (A) with redshift, except per- 
haps for a slight tendency for the photometric redshifts 
to be underestimated at z > 1.3. The rms values are all 
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based on a total of 433 galaxies in K20 (270) and F0RS2 
(163), with the object at z = 2.8 excluded (see 4.2b). 

For redshifts estimated with Bayesian priors, we find 
that the fraction of catastrophic outliers, rj, ranges from 
0.024 to 0.10, depending on the subsamplc considered. As 
described in Benitez (2000), the value of the ODDS pa- 
rameter is a reasonable indicator of the reliability of the 
photometric redshift, galaxies with larger ODDS values 
have a smaller rate of outliers. We have measured a{A) 
using only galaxies with jAj < 0.2 (i.e., excluding the out- 
liers), and find 0.049 (F0RS2) and 0.046 (K20) from the 
BPZ method. The outlier-clipped rms is similar when us- 
ing the conventional method without priors, but the 
failure rate, 77, from this method is much larger. 

We have estimated redshifts using different combina- 
tions of the available photometric data in CDF-S, and find 
no significant differences in the cr(A) values. In particu- 
lar, wc have compared the performance with and without 
using the ACS photometry, and found no significant dif- 
ference in the overall result. The deep ACS photometry, 
however, may well be important at fainter magnitudes, be- 
yond the limits of the spectroscopic sample. 18% of the 
F0RS2 galaxies, and 13% of the K20 galaxies, are un- 
detected (< 3a) in the relatively shallow WFl U or U' 
images. However, we find no significant difference in A 
values for galaxies with and without U /U' detections down 
to the spectroscopic magnitude limit (crosses in Figure 1). 
We have also estimated photometric rcidshifts using the 
WFI and SOFI data alone, which cover a wider CDF- 
S area, including areas not imaged by ACS. We measure 
(t(A) = 0.11 (for Rab < 25). 

Figure 2 plots the redshift errors A versus galaxy mag- 
nitudes in the R and Ks bands. There is no strong trend in 
fT(A) with magnitude, except at the faintest iiTg-band lim- 
its, Ks > 22.5. The outlier fraction r] increases from 0.03 
(K,, < 20) to 0.05 (20 < Ks < 22) to 0.11 (22 < Ks < 24). 
The photometric redshifts arc also, on average, slightly 
underestimated for Kg > 22. At these faint magnitudes, 
many of the galaxies are only poorly detected, if at all, 
in the wide-field (but relatively shallow) SOFI JHKg im- 
ages, causing an increase in redshift errors. 

Results from Table 1 indicate that the photometric red- 
shift performance is correlated with the ODDS parameter, 
and is better for higher ODDS values (ODDS > 0.99). It 
is important to remember that the galaxies in the spectro- 
scopic sample are relatively bright compared to the large 
majority of objects in the GOODS fields. The ODDS pa- 
rameter, perhaps combined with other indicators, (e.g., 
the number of passbands in which the galaxy is signifi- 
cantly detected; availability of ?7-band and near-IR data; 
the width of the redshift probability distribution), offers a 
useful metric for the likely reliability of photometric red- 
shifts at magnitudes fainter than the limit of our spectro- 
scopic test sample. 

4.2. Special Galaxy Populations 

One of the main aims of the GOODS project is to iden- 
tify and study different populations of galaxies, such as 
EROs and AGNs. In this section wc carry out an analysis 
of the reliability of photometric redshifts to those; objects, 
a). EROs— The GOODS ERO sample is selected to 
have (i? — K)ab > 3.35 and is complete to Kab < 22 



mag. A total of 66 EROs have spectroscopic redshifts in 
the combined F0RS2 (36) and K20 (30) samples. Fig- 
ure 3a compares the photometric and spectroscopic red- 
shifts for EROs. There is an excellent agreement, with 
(t(A) — 0.051, and is equally good for objects classified 
as absorption and emission line systems. Furthermore, we 
find a very small fraction of outliers [rj = 3%). This perfor- 
mance is significantly better than that for the galaxy sam- 
ple as a whole. This may not be surprising: in general, red 
galaxies have stronger features in their broad band spec- 
tral energy distributions (breaks, curvature) than do blue 
ones. However, it is a very helpful result, because EROs 
are among the most difficult galaxies for spectroscopic ob- 
servations. 

The ERO population is known to consist of high rod- 
shift {z ~ 1) ellipticals and dusty starbursts (Cimatti et 
al 2002c). However, the starburst galaxy templates from 
Kinney et al (1996) that are used for photometric redshift 
estimation are not significantly reddened, and certainly do 
not match the colors of EROs. Photometric redshifts for 
the majority of the EROs, therefore, are derived from the 
elliptical galaxy spectral template, regardless of the true 
nature of the galaxies. Moustakas et al. (2003) find that 
~ 40% of the GOODS EROs are morphologically early- 
type galaxies; the rest are either disk galaxies or irregular 
systems. However, they also find that the broad band 
SEDs of the EROs in different morphological subclasses 
are virtually indistinguishable. Thus, wo expect compara- 
bly good photometric redshift estimates for most EROs, 
regardless of their intrinsic nature. 

b). AGNs — 48 galaxies in the spectroscopic K20 (31) 
and F0RS2 (17) samples have X-ray detection (Alexander 
et al 2003). Figure 3b compares photometric and spectro- 
scopic redshifts for these X-ray sources. The rms scatter, 
excluding the object at ^; = 2.8, which is confirmed to 
be a QSO; is (t(A) = 0.104, with an oulier fraction of 
77 = 0.11. Excluding the five outliers reduces the rms scat- 
ter to (t(A) = 0.042. The majority of the X-ray sources at 
higher redshifts here are AGNs although some are likely 
to be X-ray starbursts. Six of these sources are spectro- 
scopically confirmed as AGNs (crosses in Figure 3b). The 
scatter and the outlier fraction for the X-ray sources (i.e. 
AGNs), is similar to that for normal galaxies in Figure 1. 
We tried using an independent set of observed AGN SEDs 
as templates, but did not significantly improve cr(A). Us- 
ing the conventional method without priors increases 
fT(A) value to 0.43. 

5. CONCLUSIONS 

Combining multi-waveband ground-based and ACS pho- 
tometric data (with up to 18 photometric passbands per 
object), we have estimated photometric redshifts for 433 
galaxies with spectroscopic redshifts in the GOODS CDF- 
S field, using the photometric redshift method of Benftez 
2000 which incorporates Bay(«ian niagnitude/redshift pri- 
ors. We find an overall performance (t(A) = 0.11, with 
better performance for some subsamples. The fraction 
of catastrophic redshift outliers is <10%, and is substan- 
tially smaller for galaxies with high values of the Bayesian 
ODDS parameter. We see no strong trend in the per- 
formance of the photometric redshifts versus magnitudes, 
except for an increase in the outlier rate. Employing a red- 
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Table 1 

Comparison between photometric and spectroscopic redshifts 

F0RS2 K20 
Riim ODDS £7(A) 7] (t(A) t] n cr(A) 77 cr(A) rj n 
BPZ BPZ 



26.00 


all 


0.107 


0.098 


0.220 


0.11 


163 


0.072 


0.048 


0.140 


0.059 


270 


25.00 


all 


0.114 


0.099 


0.110 


0.110 


141 


0.072 


0.046 


0.142 


0.058 


261 


24.50 


all 


0.111 


0.100 


0.113 


0.100 


111 


0.073 


0.049 


0.146 


0.049 


246 


25.00 


= 1.00 


0.082 


0.038 


0.083 


0.038 


78 


0.065 


0.024 


0.068 


0.030 


165 


25.00 


> 0.99 


0.095 


0.062 


0.098 


0.062 


113 


0.064 


0.028 


0.067 


0.028 


246 


25.00 


> 0.95 


0.100 


0.073 


0.102 


0.073 


124 


0.070 


0.039 


0.073 


0.043 


255 



shift/magnitude prior in this process seems to be crucial 
in reducing the scatter between photometric and spectro- 
scopic redshifts. 

We have applied the method to two sub-samples of 
galaxies of particular interest: EROs and AGNs. The re- 



sults for EROs are more accurate ((t(A) = 0.051) than for 
faint galaxies in general. They are somewhat less accurate 
for the X-ray sources (i.e. AGNs)- ((j(A) = 0.104), but 
good enough to be useful for many applications. 
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Fig. 1. — Comparison between photometric and spectroscopic redshifts for the K20 (open circles) and FORS (filled circles) samples. Galaxies 
undetected (very faint) in U and U' bands in both samples are plotted by crosses. 
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Fig. 2. — Dependence of the accuracy of photometric redshifts to brightness of galaxies in Rab and Ksj^g- Objects from both the K20 
(open circles) and FORS (filled circles) samples are included. 
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Fig. 3. — Compajrison between photometric and spectroscopic redshifts for (3a). EROs from K20 (open circles) and FORS (filled circles). 
(3b). X-ray sources (i.e. AGNs) from K20 (open circles) and FORS2 (filled circles), with spectroscopically confirmed AGNs (crosses). 



